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ABSTRACT 

Approximate  analytic  solutions  have  been  obtained  to  the  lateral 
motio.n  of  re-entry  vehicles  during  constant  altitude  glide.  The  con¬ 
stant  altitude  glide  is  maintained  by  banking  the  vehicle  about  its 
velocity  vector  at  a  fixed  angle  of  attack.  Such  a  glide  mode  is  often 
used  in  design  studies  to  achieve  lateral  range  but  no  simple  analytic 
solutions  have  been  available  up  to  now.  The  approach  taken  here  is 
to  first  treat  the  "flat  earth"  case  and  then  include  the  spherical  earth 
correction.  Simple  expressions  have  been  obtained  for  the  turn  angle, 
lateral  range  angle  and  the  down  range  angle.  A  comparison  with  exact 
numerical  solutions  shows  that  the  approximate  solutions  are  sufficiently 
accurate  for  preliminary  design  purposes. 

The  approximate  solutions,  although  simple,  contain  all  the  charac¬ 
teristics  of  the  glide  mode  in  question.  The  results  show  that  the  lateral 
range,  being  a  function  of  W/(C^A  p),  exhibits  a  maximum  for  a  given 
L/D  and  initial  velocity  for  the  spherical  earth  case.  In  the  flat  earth 
case,  the  maximum  exists  for  suborbital  and  orbital  initial  velocities 
and  ceases  to  exist  above  a  certain  initial  velocity  depending  on  L/D. 

The  down  range  also  exhibits  a  maximum  for  a  given  L/D  and  initial 
velocity  provided  that  the  initial  velocity  is  suborbital.  The  spherical 
earth  correction  on  the  lateral  and  down  range  has  been  shown  to  be 
small  for  low  W/(C^A  p),  L/D  and  initial  velocities  and  becomes  signifi¬ 
cant  when  these  parameters  are  high..  All  these  findings  agree  with 
exact  numerical  solutions. 
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I.  INTRODUCTION 

In  many  cases,  a  space  vehicle  entering  the  earth's  atmosphere  requires 
lateral  range  capability  for  landing  at  a  specific  site.  To  gain  lateral  range, 
the  vehicle  must  be  steered  out  of  the  original  orbit  plane  and  as  a  result  one 
needs  an  understanding  of  the  lateral  motion  of  such  a  vehicle.  Slye  (Reference  1) 
has  presented  an  analytical  method  for  studying  the  lateral  motion  of  re-entry 
vehicles  flying  equilibrium  glide  trajectories.  While  this  class  of  trajectories 
is  of  great  interest  and  importance  in  designing  a  re-entry  flight  path,  its 
practical  use  is  limited  to  re-entry  from  orbital  speed.  When  the  vehicle  enters 
the  atmosphere  at  superorbital  speeds,  a  transition  phase  is  required  to  bridge 
the  gap  between  the  superorbital  initial  phase  and  the  suborbital  equilibrium 
glide  phase.  One  type  of  trajectory  capable  of  carrying  out  this  transition  is  the 
constant  altitude  glide.  Constant  altitude  glide  has  drawn  the  attention  of  many 
designers;  for  example,  Smith  and  Menard  (Reference  2)  used  a  constant  altitude 
glide  to  maneuver  for  range  near  orbital  speed.  Ferri  and  Ting  (Reference  3) 
also  studied  constant  altitude  glide  to  reduce  the  second  peak  deceleration  and  to 
control  the  vehicle  from  skipping.  These  authors,  however,  did  not  consider 
the  advantage  of  obtaining  lateral  range  by  banking  the  vehicle  while  maintaining 
the  constant  altitude  glide.  Banking  the  vehicle  not  only  produces  lateral  range 
but  also  ease  the  problem  nf  designing  a  thermal  protection  system  due  to  the 
fact  that  only  one  side  of  the  vehicle  is  exposed  to  the  free- stream.  This  sub¬ 
ject  has  been  discussed  by  Gervais,  duPont  and  Lowe  (Reference  3). 

While  the  constant  altitude  glide  is  of  great  practical  interest,  it  is  not 
without  problems  to  get  a  vehicle  onto  such  a  glide  from  the  initial  phase  of  the 
re-entry,  especially  when  the  speed  is  superorbital.  It  can  be  readily  appre¬ 
ciated  that  lift  modulation  in  the  initial  phase  is  required.  Wang  and  Chu 
(Reference  4)  have  attempted  to  analytically  solve  the  problem  of  re-entry 
trajectory  with  variable  lift  of  a  simple  type.  They  have  demonstrated  that  by 
varying  the  parameters  in  the  lift  program,  one  is  indeed  able  to  make  a  smooth 
transition  from  the  initial  plunge  to  a  constant  altitude  glide.  This  can  be  done 
either  by  pitch  control  or  by  bank  control.  In  the  present  paper,  it  will  be 
assumed  that  the  vehicle  is  already  at  the  beginning  of  a  constant  altitude  glide 
and  the  initial  phase  of  the  re-entry  trajectory  is  beyond  the  scope  of  this 
analysis . 
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The  particular  maneuver  to  be  dealt  with  here,  namely,  constant  altitude 
glide  by  banking  at  a  fixed  angle  of  attack,  has  been  studied  by  Skulsky  (Refer¬ 
ence  5).  He  presented  exact  numerical  solutions  to  the  equations  of  motion. 
Later,  Wang  and  Skulsky  (Reference  6)  discussed  the  characteristics  of  the 
lateral  range  obtained  during  such  a  glide.  Their  discussion  was  based  on  the 
numerical  solutions  presented  in  Reference  5.  The  present  paper  presents, 
however,  analytic  solutions  to  the  problem  using  a  basic  approximation.  The 
approximation  appears  to  be  crude  at  first  but  the  results  turn  out  to  be  quite 
satisfactory.  The  “flat  earth1*  solution  will  be  given  in  the  first  part  of  the 
paper  and  in  the  second  part  the  spherical  earth  case  will  be  dealt  with.  A 
summary  of  the  approximate  solutions  is  given  in  Appendix  A. 
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II.  EQUATIONS  OF  MOTION 


The  three-dimensional  equations  of  motion  of  a  hypervelocity  vehicle  flying 
over  a  spherical  earth  have  been  derived  by  London  (Reference  7),  by  Love  and 
Neustadt  (Reference  8),  and  by  Wang  (Reference  9).  In  Reference  9,  the  lateral 
centrifugal  force  term  has  been  discussed  extensively  with  regard  to  its  origin 
and  physical  identity.  It  is  not  considered  necessary  to  include  the  derivation 
of  these  equations  in  this  paper.  For  constant  altitude  glide  these  equations 
reduce  to 


f  ZW  \  dV 

vcda  py dt 


.,2  'L\ 

V  cos 


0  =  ! 


2W 


cda 


11  - 


/v 

IV" 


(2) 


d\l/  _  (L\  Vg  sin  0  Vg  cos  i|/  tan  X. 

dF  '  !D'  /  2W  •  7^2 

(cDA  p; 


Here  <p  is  the  bank  angle  measured  from  the  vertical,  tJ/  is  the  turn  angle 
measured  from  the  original  heading  and  \  is  the  lateral  range  angle  defined 
as  the  ratio  of  the  lateral  range  to  the  radius  of  the  earth  (equal  to  3440 
nautical  miles  in  the  calculations  shown  later).  These  angles  are  shown  in 
Figure  I.  The  other  symbols  have  their  usual  meaning.  Together  with  these 
three  equations  of  motion,  there  are  two  kinematic  equations: 


d\  Vg  sin  i); 

dF  '  ^2 


S 


(4) 


- 

dt  " 


Vg  cos  ij/ 

V2  cos  \ 
s 


(5) 


Here  p  is  the  downrange  angle  as  shown  in  Figure  1.  It  is  to  be  noted  that 
Equations  (1),  (2),  and  (3)  are  applicable  for  a  point  mass.  A  point  mass 
vehicle  is  assumed  in  this  analysis.  The  orbital  speed  Vg  is  assumed 
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Figure  la.  Definition  of  Bank  Angle  0. 

(Velocity  vector  pointing  into  the  paper) 


Figure  lb.  Definition  of  Turn  Angle  ty,  Lateral  Range  Angle 
and  Down  Range  Angle  p. 
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constant  and  equal  to  26,  000  ft/ sec  in  the  numerical  calculations  given  later 
in  the  paper. 

If  we  write 


B 


W 

cda  p 


(6) 


and  transform  the  independent  variable  from  time  t  to  velocity  V  through 
Equation  (1),  we  obtain 

-  V  sin  0  -  cos  ljj  tan  X  (7) 

^  s 


-  V 


dX 

HV 


(8) 


y.  dji  _  2B  cos  ljj 
V  3V  "  cos  X 


Denoting  In  V  =  X,  these  equations  become 


dt|r 

3X  = 


L  .  ^  2B  ,  ,  x 

sm  0 - j  cos  f  tan  X 


(9) 


(10) 


dX 

dX 


2B  .  . 

— sin  i[r 


(ii) 


dp  2B  cos  ljj  ,  ^ 

dX  ~  yZ  cos  X  ' 

s 

Equations  (10),  (11),  and  (12)  determine  the  unknowns  tlr,  X,  and  u  in  terms 
of  V.  The  bank  angle  0  can  be  readily  written  in  terms  of  V  by  Equation  (2) 
and  the  velocity  V  can  be  found  as  a  function  of  t  by  integrating  Equation  (1). 
The  main  problem  here  is  to  solve  the  three  nonlinear  differential  Equations 
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(10),  (11),  and  (12).  One  can  readily  appreciate  the  difficulties  associated 

with  solving  these  equations  analytically.  The  approach  taken  here  is  to  treat 

the  "flat  earth"  case  first  and  then  take  up  the  spherical  earth  case.  Flat 

earth  means  a  condition  when  the  second  term  on  the  right  hand  side  of 

Equation  (10)  can  be  neglected.  This  term  is  proportional  to  the  lateral 

centrifugal  force.  If  the  earth  is  truly  flat,  this  term  will  of  course  drop  out 

2 

but  so  will  the  (V/Vg)  term  in  Equation  (2).  In  the  flat  earth  solution  pre¬ 
sented  below,,  Equation  (2)  remains  the  way  it  is  and  only  Equation  (10)  is 
simplified  by  neglecting  the  second  term  on  the  right  hand  side.  By  so  doing, 
it  does  not  result  in  a  true  flat  earth  as  the  name  implies. 


Page  7 


III.  FLAT  EARTH  SOLUTION 

For  the  flat  earth  case,  although  the  equations  to  be  solved  are  much 
simplified,  they  are  still  difficult  to  handle  analytically  because  of  the  trigono¬ 
metric  functions  involved.  One  certainly  can  eliminate  sin  0  by  Equation  (2) 
and  integrate  Equation  (10)  and  in  principle  it  is  integrable  for  the  case  of  flat 
earth.  But  the  resultant  expression  for  will  be  very  complicated.  In  fact, 
this  has  been  done  by  Skulsky  (Reference  10),  and  one  still  has  to  use  high 
speed  computing  machines  to  get  a  numerical  value.  As  for  X  and  p,  numeri¬ 
cal  integration  is  unavoidable.  What  is  done  here  is  to  further  simplify  Equation 
(10)  by  making  the  following  approximation: 


sin  0=1  (13) 

Let  us  investigate  this  approximation.  First,  rewrite  Equation  (2)  as  follows: 


,[ 

7\m2  ' 

sin  0  =  “J 

i  -c2 

V ~V7  "  \ 

(14) 


where 


C 


2W/(ClA  p  V2) 


(15) 


Equation  (14)  is  plotted  in  Figure  2.  It  is  seen  in  Figure  2  that  for  small  values 

of  C,  sin  0  is  practically  equal  to  unity  for  a  wide  range  of  V/Vg.  For  example, 

when  the  initial  velocity  is  V\/Vg  =  1-2  and  C  =  0.  1,  sin  0  is  practically  unity 

down  to  V/Vg  =0.5  and  rapidly  decreases  to  zero  at  V\ /Vg  =  0.  3.  Integrating 

Equation  (10)  (when  the  second  term  on  the  right  hand  side  is  omitted  for  flat 

earth)  is  effectively  to  find  the  area  under  each  sin  0  curve  for  various  C. 

The  approximation  sin  0=1  will  give  the  area  of  the  rectangle  bounded  by 

V./V  and  V./V  where  V-  is  the  final  velocity  that  is  at  the  end  of  the  constant 
l  s  f  s  f  7 

altitude  glide  where  0=0.  This  area  will  not  be  much  larger  than  that  under 
the  actual  sin  0  curve.  For  high  C,  the  error  should  increase.  For  low  initial 
velocities,  the  error  should  also  increase.  For  example,  when  W/Vg  =  0.77 
and  C  =  0.8,  there  is  no  point  where  sin  0  is  close  to  unity.  Therefore, 
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the  approximation  is  the  worst  at  high  C  and  low  initial  velocities.  The 
numerical  examples  given  later  in  the  paper  bear  this  out  but  the  inaccuracy 
in  \  and  \i  is  surprisingly  low. 

By  this  approximation  sin  0=1,  Equation  (10)  for  the  flat  earth  case 
reduces  to 


d*o 

■ax 


L 

D 


(16) 


Here  the  subscript  0  denotes  flat  earth  solution.  The  initial  condition  is 
l| Jq  =  0  at  .  Integrating  Equation  (15)  gives  the  following  simple  result. 

*0  =  fe<xi  -  X>  (17) 

or 

*0  =  pln(V./V)  (18) 


By  Equation  (17),  Equation  (11)  can  be  integrated  with  the  initial  condition  that 
\«  =  0  at  X.  and  the  result  is 

0  l 

\0  =  C  ( 1  -  cos  t0)  (19) 


At  the  end  of  the  constant  altitude  glide,  the  bank  angle  0  vanishes;  that  is  to 
say,  the  vehicle  has  reached  the  saturation  point  and  cannot  maintain  a  constant 
altitude  unless  by  changing  the  angle  of  attack.  Since  the  flight  mode  considered 
here  is  constant  altitude  glide  by  banking  at  a  fixed  angle  of  attack,  the  final 
velocity  must  be  determined  by  the  condition  0=0.  Therefore,  by  Equation 
( 14),  one  has 

V 

s 


(20) 
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The  final  turn  angle  and  the  final  lateral  range  angle  are  respectively 


\i/ 


Of  " 


^of  =  c  U  '  cos  v0f) 


(21) 

(22) 


At  this  point,  it  is  worthwhile  to  compare  these  simple  solutions  with 
exact  numerical  solutions.  The  exact  numerical  solutions  ar.e  taken  from 
References  5  and  10.  Reference  10  gives  the  flat  earth  solutions  and  should 
serve  as  a  check  against  the  sin  <£  =  1  approximation.  Reference  5  gives  the 
spherical  earth  solution  and  should  indicate  the  magnitude  of  the  spherical 
earth  effect. 

Figure  3  is  a  plot  of  ^^/(L/D)  as  a  function  of  W/(C^A  p)  for  initial 
velocities  32,  000,  26,  000,  and  20,  000  ft/sec.  Equation  (2  1)  predicts  too  high 
a  turn  angle  as  expected.  The  percentage  error  increases  with  increasing 
W/(CT  A  p)  for  a  given  initial  velocity  and  decreases  with  increasing  initial 

Jj 

velocity  for  a  given  W/(CT  A  p).  These  are  the  expected  nature  of  the  approxi- 

l—i 

mation.  The  absolute  magnitude  of  the  error  in  the  quantity  ^^/(L/D)  is  of 
the  order  of  2  to  about  6  degrees  in  the  range  compared.  Generally  speaking, 
the  accuracy  of  the  simple  approximate  solution  is  satisfactory;  it  correctly 
predicts  the  general  trend  and  correctly  establishes  the  functional  dependence 
of  upon  B,  L/D  and  V..  Figures  4a,  4b  and  4c  show  the  same  comparison 
plotted  in  a  slightly  different  form.  Included  in  these  figures  are  the  exact 
numerical  solutions  for  spherical  earth,  which  will  be  discussed  later  in 
Section  IV. 

A  comparison  of  the  lateral  range  is  given  in  Figures  5a,  5b  and  5c.  For 
all  three  initial  velocities,  the  approximate  solution,  Equation  (22),  agrees 
with  the  exact  numerical  solutions  for  flat  earth  well  within  10  per  cent.  In 
some  cases  when  the  initial  velocity  is  high,  the  present  approximate  solution 
is  almost  exact  and  so  is  it  for  low  W/(C^A  p).  The  general  characteristics 
of  the  lateral  range  as  a  function  of  B,  L/D  and  V.  including  the  peaks  are 
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Figure  4a.  Comparison  of  Final  Turn  Angle  (V.  =  Z0,  000  ft/sec). 
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Figure  4b.  Comparison  of  Final  Turn  Angle  (V\  -  26,  000  ft/ sec). 
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Figure  4c.  Comparison  of  Final  Turn  Angle  (V\  =  32,  000  ft/sec). 
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Figure  4c.  Comparison  of  Final  Turn  Angle  (V\  =  32,  000  ft/sec). 
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of  Lateral  Range  (V^ 


1000  b - - - -L  ,  2.0  - - -  -  EQUATION  (22) 
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Figure  5b.  Comparison  of  Lateral  Range  (V.  =  26,  000  ft/ sec 


Comparison  of  Lateral  Range  (V.  =  32,  000  ft/sec) 
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well  preserved  in  the  approximate  solution.  The  striking  feature  that  one  finds 
in  these  comparisons  is  that  the  accuracies  of  predicting  ^  and  are  not 

at  all  the  same.  The  approximate  solution  predicts  the  lateral  range  with  a 
much  higher  accuracy  than  it  does  the  turn  angle.  This  can  probably  be 
explained  in  the  following  manner.  For  a  given  initial  velocity  and  W/(C^A  p), 
the  exact  flight  path  might  be  represented  by  the  solid  curve  of  Figure  6  while 
the  approximate  solution  is  indicated  by  the  dashed  line.  Point  A  might  repre¬ 
sent  the  end  of  the  turn.  The  lateral  range  in  both  computations  might  be  very 
close  to  each  other  while  the  approximate  solution  predicts  a  much  higher  turn 
angle  than  the  exact  numerical  solution.  This  explanation  is  fully  supported  by 
the  downrange  prediction.  The  approximate  solution  always  predicts  a  shorter 
downrange  as  will  be  seen  later. 

Let  us  now  return  to  the  development  of  the  approximate  solution.  Equations 
(2  1)  and  (22)  can  be  used  to  determine  the  maximum  lateral  range  and  the  con¬ 
dition  under  which  it  occurs.  Differentiating  Equation  (22)  with  respect  to  C 
and  setting  it  equal  to  zero,  one  gets 


d\ 


Of 


dill 


~sc~ 


=  (1  -  cos  4q^)  +  C  sin  y 


Of 


'Of  ~dcT 


r  0 


where  d^^/dC  can  be  found  from  Equation  (21)  as 


dijr 


Of 


1  L 


1 


dC  2  D  C  (1  +  C) 

Combining  Equations  (24)  and  (23),  one  gets 


1  -  cos  ill 


1  L  Sin  ^  Of 


Of 


7  D 


1  +  C 


(23) 


(24) 


(25) 


Here  the  #  quantities  denote  the  condition  where  the  maximum  lateral  range 
occurs.  Using  the  trigonometric  relation  of  half  angles  and  eliminating  i|i  ^ 
by  Equation  (21),  Equation  (25)  finally  reduces  to 

T - ; — 

Dr.  Frank  Billet  analytically  investigated  the  errors  and  showed  similar 
results. 


LATERAL  RANGE 
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■  1  '  REPRESENTING  EXACT  FLIGHT  PATH 

—  —  —  REPRESENTING  APPROXIMATE  FLIGHT  PATH 


Page  20 


V. 


1  L\  ,  '  wi  L  .  1 

1  Di  ln  ,V  ’V1  +  C* 

s  v 


arc  tan 


1  h 

2  D 

1  T~C*! 


(26) 


The  maximum  lateral  range  can  be  readily  found  to  be 


r.v.  * 


e  •  ,  L  T  i  ;  T  i  s  !  ,  1 

sin  ^  In  .j  y-;  -1  +  £*./• 
/ 1  Li .  s 


l0f  “  1VI  D/ 


1  + 


C* 


(27) 


These  equations  have  been  used  to  compute  the  maximum  lateral  range  and 
where  it  occurs.  The  results  are  listed  in  Table  1.  Included  in  Table  1  are 
also  the  exact  numerical  solutions.  The  comparison  is  extremely  favorable. 

It  is  found  that  for  a  given  L/D  there  exists  an  initial  velocity  beyond  which 
there  is  no  maximum  for  a  flat  earth.  This  is  illustrated  in  Figure  7  in  which 
Equation  (26)  is  plotted. 

The  down  range  angle  can  now  be  integrated.  Dividing  Equation  (12)  by 
Equation  (11)  one  gets  for  flat  earth, 


dp 

dX 


0 

0 


cot 


co  s 


(28) 


By  Equation  (19),  one  can  write 

C  -V 

cot  *  =  _ ^=r  (29) 

'Jc2  -  (C  -  \Q)2 

1  2 

The  1  / cos  term  can  be  approximated  by  l  +  ^  \q.  This  approximation  is 
accurate  within  5  per  cent  up  to  =  40  degrees  which  means  lateral  range  = 
2400  nautical  miles.  In  Figures  5a,  5b  and  5c,  the  lateral  range  is  smaller 
than  2400  nautical  miles  for  all  the  conditions  presented.  If  the  lateral  range 


Table  1.  Maximum  Lateral  Range  and  B*  (= 
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Figure  7.  Comparison  of  C*,  "Flat  Earth 
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is  of  the  order  of  1000  nautical  miles,  the  approximation  has  an  inaccuracy  of 
about  1  per  cent.  Together  with  this  approximation,  Equations  (28)  and  (29) 
combine  to  yield 


<^0 

W~0 


<c  - 

-\Jc2  -  (C  -  \Q)2 


(30) 


The  initial  condition  is  =  0  at  \q  =  0.  One  finally  gets  by  integrating 
Equation  (30) 


.  c2  .  cvofL  g3  c3  I  .  /.  ''of y  n 

^of  =  V1  +  ~r +  2  /  G  "  ~r  +  ~r  iarc  sin  l1  ■  ~c~) - 1 


(31) 


where 


G  =  -/C2  -  (C  -  \Qf)2 


(32) 


In  terms  of  q£*  Equation  (31)  becomes 

3  f  1  1  s^n  tgfl 

^0f  =  C  sin  ^ Of  +  C  sin  *0£  -  7  *0f  ■  T  sin  2  *0f - 6 - 


(33) 


For  practical  purposes,  Equation  (33)  can  be  approximated  by 


M-of  =  C  sin  i) 


Of 


(34) 


Equation  (34)  agrees  with  Equation  (33)  within  1  per  cent  for  lateral  range  up 
to  1000  nautical  miles  and  within  3. 5  per  cent  up  to  2000  nautical  miles.  Note 
that  Equation  (34)  would  be  the  result  of  the  integration  if  cos  =  1  were 
assumed.  The  downrange  is  plotted  in  Figures  8a,  8b  and  8c.  Comparison 
with  the  exact  numerical  solutions  show  that  the  approximate  solution  predicts 
shorter  downrange  as  depicted  by  the  schematics  of  Figure  6,  but  the  accuracy 
is  very  high.  In  these  plots,  Equation  (34)  is  used  except  for  a  few  points  where 
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Figure  8a.  Comparison  of  Down  Range  (V.  =  20,  000  ft/sec). 
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Figure  8c.  Comparison 
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the  lateral  range  is  above  1500  nautical  miles.  It  is  observed  that  for  a  sub¬ 
orbital  initial  velocity,  there  exists  a  maximum  down  range.  No  maximum 
exists  for  orbital  or  superorbital  initial  velocities  because  the  vehicle,  at 
these  speeds,  can  always  fly  at  a  sufficiently  high  altitude  and  go  around  the 
earth  as  many  times  as  desired.  This  would  produce  no  maximum  on  down- 
range.  It  is  also  noted  that  the  maximum  down  range  and  maximum  lateral 
range  do  not  occur  at  the  same  altitude  for  a  given  vehicle.  The  maximum 
down  range  occurs  at  a  higher  altitude.  One  can  readily  find  the  maximum 
down  range  at  suborbital  initial  velocities  by  differentiating  Equation  (34)  and 
following  the  usual  procedures.  The  result  is  given  in  Appendix  A. 

For  completeness,  the  total  flight  time  is  given  below.  This  is  obtained 
by  integrating  Equation  (1)  with  the  initial  condition  t  =  0  at  V.  and  the  final 
condition  t  =  t^  at  V^. 


This  essentially  completes  the  development  of  the  flat  earth  solution.  It 
must  be  pointed  out,  however,  that  the  solutions  presented  above  must  be  used 
with  caution  when  W/(C^A  p)  becomes  very  small.  In  the  previous  figures, 
the  lateral  and  down  range  have  been  extended  to  zero  at  W/(C^A  p)  =  0  by 
smooth  curves.  This  was  done  as  a  matter  of  convenience.  Strictly  speaking, 
this  is  incorrect.  As  one  can  readily  see  that  for  a  given  L/D  and  initial 
velocity,  the  down  range  and  lateral  range  expressions  Equations  (22)  and  (34) 
can  be  combined  to  yield  an  equation  of  a  circle  with  its  center  at  X  =  C  and 
radius  equal  to  C  .  When  C  is  extremely  small,  although  this  may  never 
become  the  case  in  practice,  the  vehicle  can  circle  around  the  X  =  C  location 
without  gaining  much  lateral  range  and  down  range.  The  lateral  range  and 
down  range  curves  must  have  an  oscillatory  nature  at  very  small  W/(C^A  p)  . 
A  more  detailed  discussion  is  given  in  Appendix  B. 
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IV.  SPHERICAL  EARTH  SOLUTION  - 
SPHERICAL  EARTH  CORRECTION 

As  can  be  seen  from  Figures  4,  5  and  8,  the  exact  numerical  solutions  for 

spherical  earth  deviate  only  slightly  from  those  for  flat  earth  when  the  initial 

velocity  is  low.  For  example,  when  W  =  20,  000  ft/sec,  one  can  hardly  see  any 

difference.  This  is  true  for  the  turn  angle,  the  lateral  range  as  well  as  the 

down  range.  When  =  26,  000  ft/ sec,  the  spherical  earth  solution  begins  to 

show  some  difference.  The  difference  on  lateral  range  is  about  10  per  cent  at 

6  2  2 

W/[CdA  p)  =  1000  x  10  ft  /sec  and  L/D  =  2.  The  spherical  earth  solution 

shows  higher  difference  at  higher  initial  velocities.  For  example,  when  W  = 

32,  000  ft/sec,  the  spherical  earth  solutions  deviate  from  the  flat  earth  solutions 

6  2  2 

on  lateral  range  by  about  20  per  cent  at  W^C^A  p)  =  1000  x  10  ft  /sec  and 
L/D  =  2.  This  order  of  magnitude  of  the  error  produced  by  the  flat  earth  solu¬ 
tion  is  too  high  to  be  ignored  and  must  be  taken  into  consideration.  On  the  other 
hand,  an  error  of  20  per  cent  or  so  can  still  be  treated  mathematically  as  a 
second  order.  Thus,  the  second  term  on  the  right  hand  side  of  Equation  (10) 
can  be  treated  as  a  small  perturbation  on  the  first  term.  Accepting  this  idea, 
one  is  justified  to  speak  of  a  spherical  earth  correction  and  treat  it  accordingly. 

In  the  case  of  an  equilibrium  glide  trajectory,  Slye's  solutions  (Reference  1) 
are  for  flat  earth.  The  spherical  earth  correction  in  that  case  has  been  dis¬ 
cussed  by  Jackson  (Reference  11)  and  by  London  (Reference  12).  London 
advanced  an  argument  that  the  proper  criterion  for  determining  the  validity  of 
the  flat  earth  solution  is  the  relative  magnitude  of  the  integrated  lateral 
centrifugal  force  to  the  integrated  lateral  aerodynamic  force,  while  the  magni¬ 
tude  of  the  lateral  range  is  not  in  itself  a  sufficient  criterion.  In  the  present 
case,  we  also  observe  that  it  is  not  sufficient  to  say  when  \  is  small,  the 
lateral  centrifugal  force  term  can  be  dropped  from  Equation  (10).  If  one  tries 

to  do  this,  he  will  run  into  difficulties,  because  the  first  term  (L/D)  sin  (p  tends 

2 

to  zero  toward  the  end  of  the  glide  while  the  second  term  (2B  cos  f  tan  \)/Vg 

can  be  very  large.  However,  if  one  compares  the  integrated  values  of  the  term 

2 

(L/D)  sin  <f>  and  the  term  (2B  cos  i|/  tan  \)/Vg  ,  he  will  find  that  for  most  cases 
the  latter  will  be  smaller  than  the  former.  If  the  latter  is  so  small  that  one  can 
neglect  it,  the  flat  earth  solution  becomes  valid. 
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The  spherical  earth  correction  for  constant  altitude  glide  has  been  dis¬ 
cussed  in  Reference  6  where  an  empirical  correlation  was  given  for  L/D  =  2. 
What  will  be  done  here  is  to  find  an  analytic  expression  for  this  correction 
v/ithin  the  present  framework  of  approximation.  Let  us  write 

♦  =  *0  +  (36) 

thus 


d*  d*0Ad*l 

HX  "  + "XT 


Combining  Equations  (10)  and  (37),  one  has 


dt|f 


o  d*  1 


-  *  L  .  ,  2B  ,  . 

'  "ax'  ‘  -  D  Sln  *  ~  ^7  COS  *  tan  X 

v  s 


(37) 


(38) 


By  the  sin  <p  -  1  approximation,  the  first  term  on  both  sides  of  the  equation 
cancels  each  other  and  one  gets 


di|f  i 

sx 


cos  tjf  tan  X 


(39) 


With  fairly  good  accuracy,  one  can  replace  cos  t|/  and  tan  X  by  cos  \jj q  and 
Xq,  respectively,  and  get 


ddf  j 


ZB  X, 


7T  cos  *0 


(40) 


or 


dtlt. 


C  (1  -  cos  i)(q)  cos  di|/ 


'0 


(41) 


Page  30 


The  initial  condition  is  yjr ^  =  0  and  i|/q  =  0.  Integrating  Equation  (41)  yields 

♦  if  =  -  C2  (sin  *0f  ■  7  #0f  '  T  sin  2  W  (42) 

Equation  (42)  relates  the  spherical  earth  correction  and  the  flat  earth  solution. 

2 

It  is  noted  that  il/  ^  is  proportional  to  C  and  it  is  always  negative.  Equation 
(42)  is  plotted  in  Figure  9  and  also  tabulated  in  Table  2.  Applying  this  correc¬ 
tion,  the  approximate  spherical  earth  solution  becomes 

2  11 

tf  =  *0f  -  C  (sin  *of  -  J  f0£  -  ^  Sin  2  i|r Qf )  (43) 

where  has  been  given  by  Equation  (21).  For  comparison,  Equation  (43)  is 

plotted  in  Figures  4a,  4b  and  4c  as  the  dashed  lines.  The  comparison  is  quite 
favorable . 

The  spherical  earth  correction  on  lateral  range  is  found  as  follows.  With 
Equation  (36)  one  can  write  Equation  (11)  as 


Since 


Thus, 


d X  2B  .  f ,  \ 

*  ^  **0  +  +lJ 

S 


d\  2B  ,  .  .  .  .  ,  •  .  v 

-  g~  =  — 2  (sm  ilr  0  cos  it  ^  +  cos  sm  ^) 

^  s 


is  small,  one  is  justified  to  assume 

cos  if  j  =  1  and  sin  if  ^  =  ij;  ^ 


(44) 


(44) 


(45) 


d\  _  2B 

"3X  =  7 


s 


sin  i|/r 


2B 

7 

s 


COS  llfQ 


(46) 
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Table  2.  Spherical  Earth  Correction  on  Turn  Angle. 


*0f 

♦  of /cZ 

(deg)  . 

(deg) 

0 

0 

10 

-  0.063 

20 

-  0.384 

30 

-  1.245 

40 

-  2.728 

50 

-  4.710 

60 

-  7.  175 

70 

-  9.  65 

80 

-  11.48 

90 

-  12. 33 
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Table  2.  Spherical  Earth  Correction  on  Turn  Angle. 


♦of 

*o  f/c2 

(deg) 

(deg) 

0 

0 

10 

-  0.063 

20 

-  0.384 

30 

-  1.245 

40 

-  2.728 

50 

-  4.710 

60 

-  7. 175 

70 

-  9. 65 

80 

-  1 1 . 48 

90 

-  12. 33 
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Following  Equation  (36),  one  can  write 


X  =  \q  +  X^ 


and  Equation  (46)  is  reduced  to 


dX  q  d\  j  2  3  2  B 

ST 'ST  =  ^Z  sin  *0  +  ^Z  *i  cos  ♦o 

s  s 


Cancelling  the  flat  earth  terms,  i.  e., 


dX 


0  2B  . 


■ax 


sin  i|rf 


one  gets  for  the  spherical  earth  correction 


dX  ^  2  g 

■air  =  "  ^z  *i  cos  +o 

s 


(47) 


(48) 


(49) 


(50) 


or 


3  /  *o  sin  2  ^o\ 

d\1  =  C  cos  sin  *0  +  -2-  + - 5 - ^  dt0  (51) 

The  initial  condition  is  X^  =  0  at  ^  =  0.  One  finally  gets  by  integrating 
Equation  (5 1 ) 

Vlf  =  ‘  °3  (|+^cos3  *0f  '  Icos2  *0f  ‘  Zcos  +0f  '  7* 0f8in*0f)  (52) 

Equation  (52)  is  plotted  in  Figure  10  and  also  tabulated  in  Table  3.  Note  that 

3 

X^,  being  negative,  is  proportional  to  C  .  Applying  this  correction,  the 
approximate  spherical  earth  solution  becomes 


90' 
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Figure  10.  Spherical  Earth  Correction  on  Lateral  Range  Angle. 
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Table  3.  Spherical  Earth  Correction  on  Lateral  Range  Angle 
and  Lateral  Range. 


*0f 

\1(/c3 

3 

Lateral  Range/C 

(deg) 

(radians ) 

(naut  mi) 

0 

0 

0 

10 

-  0.  00075 

-  .2.58 

20 

-  0.  00170 

5.  85 

30 

-  0.  00325 

-  11.18 

40 

-  0.  00730 

-  25.  1 

50 

-  0.  0157- 

-  54.  0 

60 

-  0.  0266 

-  91.5 

70 

-  0.0365 

-  125.5 

80 

-  0.  0458 

-  157.5 

90 

-  0.  0483 

-  166.  0 
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Equation  (53)  is  plotted  in  Figures  5a,  5b  and  5c  as  the  dashed  lines  and  is 
compared  with  the  exact  numerical  solutions  from  Reference  5.  The  comparison 
shows  that  the  corrected  \  expression  can  predict  the  lateral  range  very  closely 
for  most  cases.  One  interesting  thing  is  that  whenever  the  spherical  earth  cor- 
rection  is  negligible  in  the  exact  numerical  solutions,  it  is  also  negligible  in  the 
approximate  solutions.  When  the  spherical  earth  correction  is  large,  the 
approximate  solution  deviates  from  the  exact  numerical  solution  no  more  than 
10  per  cent  in  the  range  compared.  As  for  the  maximum  lateral  range,  it  has 
been  noted  in  the  flat  earth  solutions  that  there  is  no  maximum  beyond  a  certain 
V\.  With  the  spherical  earth  correction,  the  maximum  lateral  range  at  low  V\ 
is  not  affected  but  a  maximum  exists  at  high  V.  due  to  the  correction.  This  fact 
has  been  observed  in  the  exact  numerical  solutions.  It  can  now  be  concluded  that 
the  maxima  at  low  initial  velocities  are  caused  by  the  integrated  lateral  aero¬ 
dynamic  force  while  those  at  high  initial  velocities  are  caused  by  the  spherical 
earth  effect.  This  indicates  the  importance  of  the  spherical  earth  correction  at 
high  initial  velocities. 

The  spherical  earth  correction  on  the  down  range  angle  jjl  is  found  as 
follows.  Taking  Equation  (1Z)  and  substituting  i|r  by  i|Tq  ♦  ^  and  \  by  \ ^ , 

one  gets 

dfi  2BCO*<*0+*1) 

‘  HX  “  cos  (\Q  +  Kl)  '54) 

S 


or 


d|x  2B  608  *0 

HX  =  COS~\q 

v  s 


cos  i|r  j 
cos  \ 


-  tan  i|/q 

-  tan  Vq 


sin  i|f  ' 

sin 


(55) 


Since  and  are  small,  one  can  write 

cos  t|fj  =  cos  =  1,  sin  ^  sin  (56) 
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and  Equation  (55)  becomes 


dp  _  2B  COS  W1  -  +1  tan  ,,r0 
3X  cos  Xq  IT  -  Xj  tan  XQ 


(57) 


Since  X^  is  always  small  and  Xq  very  rarely  reaches  beyond  45  degrees,  one 
can  assume 


thus 


Xj  tan  Xq  «  1 


(58) 


2BCOS*0n  , 

T^(1  -*1 

S 


tan  t|»0) 


Writing  jjl  =  ^Lq  +  and  substituting  in  the  above,  one  gets 


d^l  _  2B  (h  Sin  M 
"  ^7  V  cos  ) 
s 

and  by  Equation  (16), 


sin  t lr0 


d^i  =  -c*i^75^d*o 


(59) 


(60 


(61) 


Using  the  approximation  cos  \q  =  1  which  has  been  proven  adequate  in  com¬ 
puting  the  down  range,  one  gets 


dpj  =  -  C  (fj  sin  tif 0  dt|/Q 


(62) 


or 


dp.  =  C  sin  ih 


'”0 

T 


sin  2 


sin  dfQ 


(63) 
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The  initial  condition  is  jjl^  =  0  at  \|/q  =  0.  Integrating  Equation  (63),  one  gets 


'If 


7" 


+  cos  Ilf 0f )  (t0£  -  sin  ||I 


Of 


\  1  *  o  | 

)  -  -g  sin  y 


Of 


(64) 


Equation  (64)  is  plotted  in  Figure  11  and  also  tabulated  in  Table  4.  It  is  noted 

3 

that  jjl  ^  is  positive  and  proportional  to  C  .  Using  this  correction,  the  approxi¬ 
mate  spherical  earth  solution  for  down  range  is 


=  C  sin  til  Qf  +  C 


1  (1  +  cos  I|r0f)  (tof  -  Sin  i|rof)  -  sin-3  * 


Of 


(65) 


The  corrected  down  range  is  plotted  in  Figures  8a,  8b  and  8c  as  the  dashed 
lines.  The  comparison  with  exact  numerical  solutions  is  quite  satisfactory. 
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Figure  11.  Spherical  Earth  Correction  on  Down  Range  Angle. 
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Table  4.  Spherical  Earth  Correction  on  Down  Range  Angle 
and  Down  Range. 


*0f 

(deg) 

“i i/c 3 
(radians) 

Down  Range/C 
(naut  mi) 

0 

0 

0 

10 

0. 00001 

0. 0344 

20 

0. 00022 

0.  7560 

30 

0. 00119 

4.  10 

40 

0. 0045 

15.50 

50 

0. 0126 

43.40 

60 

0. 0276 

95. 00 

70 

0.  0507 

174. 20 

80 

0.  0823 

283. 00 

90 

0.  1187 

407.50 
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V.  CONCLUDING  REMARKS 

An  analytic  study  has  been  made  of  the  lateral  motion  of  a  re-entry  vehicle 
gliding  at  a  constant  altitude.  The  constant  altitude  is  maintained  by  banking  the 
vehicle  about  the  velocity  vector  at  a  constant  angle  of  attack.  The  analytic 
approximate  solutions  have  been  compared  with  numerical  exact  solutions.  It 
has  been  shown  that  the  approximate  solutions  predict  higher  turn  angle,  higher 
lateral  range  but  lower  down  range.  The  accuracy  of  the  prediction  is  generally 
satisfactory.  Among  the  three  variables,  the  lateral  and  down  range  can  be  pre¬ 
dicted  well  within  10  per  cent  while  the  turn  angle,  with  lower  accuracy. 

There  are  two  important  features  about  these  approximate  solutions.  First, 
they  are  simple  and  can  be  readily  applied  to  preliminary  design  work.  Secondly, 
they  preserve  the  characteristics  and  provide  a  good  understanding  of  a  constant 
altitude  glide. 

The  results  of  the  approximate  solutions  show  (and  the  exact  numerical 
solutions  agree)  that  the  lateral  range,  being  a  function  of  W/(C^A  p)  at  a 
given  L/D  and  V\ ,  exhibits  a  maximum.  The  maximum  exists  for  all  initial 
velocities  in  the  case  of  a  spherical  earth.  For  a  flat  earth  it  exists  at  low 
initial  velocities  only.  The  spherical  earth  correction  has  been  shown  to  be 
small  for  low  W/(Cj^A  p),  L/D  and  and  become  significant  only  when  these 
parameters  assume  high  values. 


APPENDIX  A 


SUMMARY  OF  APPROXIMATE  SOLUTIONS 


I.  TOTAL  FLIGHT  TIME  (EXACT) 


LC  I  'si 

'<  '  If*'  '~ii 


where 


C  = 


2  W 


CL  AP  Vs 


II.  FINAL  VELOCITY  (EXACT) 


V,  = 


III.  FLAT  EARTH  SOLUTIONS 
Basic  approximation: 


Final  turn  angle: 


sin  0  =  1 


^of  =  n  ln  (\r  il1  +u 


n 


Final  lateral  range  angle: 


X0f  =  C  t1  '  cos  ^Of) 


Final  down  range  angle: 
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C*  for  maximum  lateral  range; 


V. 

1 

V” 

s 


arc  tan 


Maximum  lateral  range 


X 


* 

Of 


C*:  for  maximum  down  range 


1  +  c*j 


-  arc  tan 


Maximum  down  range 


►‘Of 


1  L 

Id 


IV.  SPHERICAL  EARTH  CORRECTION 


Correction  on  final  turn  angle 


♦if 


sin  2 


Correction  on  final  lateral  range 


Correction  on  final  down  range 


►‘if 


2  (1  +  cos  I|i0f)  (\(i0£  -  sin  tof)  -  F  Sin  Q{ 
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APPENDIX  B 

DISCUSSION  OF  THE  APPROXIMATE  SOLUTIONS 
FOR  VERY  SMALL  C+ 

It  has  been  mentioned  in  the  text  that  the  sin  0=1  approximation  becomes 
more  accurate  when  C  gets  smaller.  The  text  also  shows  that  the  spherical 
earth  correction  becomes  negligible  when  C  is  very  small.  Furthermore,  the 
lateral  range  angle  X  must  be  small  for  small  C  (X  is  proportional  to  C)  and 
as  a  result  the  cos  X  =  1  approximation  used  in  computing  the  down  range 
angle  becomes  more  accurate  as  C  gets  smaller.  Therefore,  for  very  small 
C  ,  Equations  (10),  (11)  and  (12)  can  be  written,  with  high  degree  of  accuracy, 
as 


ddf 

lx 

n 

(B-l) 

dX 

dX 

ii 

<  1^ 

*  t >j|W 

sin  i|r 

(B-2) 

dp 

dX 

2B 

‘  v2 

s 

COS  ilf 

(B-3) 

The  solutions  of  these  equations  are: 

*  = 

l  i 

D  ln 

(B-4) 

X  = 

C  (1  - 

COS  (||) 

(B-5) 

H-  = 

C  sin 

* 

(B  — 6) 

The  author  is  grateful  to  Dr.  Frank  Billet  for  his  helpful  suggestions 

in  this 

discussion. 
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Combining  Equations  (B-5)  and  (B-6),  one  gets 

(\  -  c)2  +  (X2  =  c2 


(B-7) 


Equation  (B-7)  represents  a  circle  with  its  center  at  \  =  C  and  its  radius 
equal  to  C  .  Such  a  circle  is  plotted  in  Figure  B-l.  The  tangent  at  any  point 
on  the  circumference  of  the  circle  makes  an  angle  \|r  with  the  X.  =  0  line.  The 
line  of  X  =  0  is  along  the  original  heading  of  the  vehicle  and  ilr  must  be  the 
turn  angle.  Figure  B-l  clearly  indicates  that  the  vehicle,  while  making  a 
constant  altitude  turn,  describes  a  circular  path.  Depending  on  the  aero¬ 
dynamics,  the  mass  characteristics  and  the  altitude,  the  vehicle  may  find 
itself  anywhere  along  the  circumference  of  the  circle  at  the  end  of  the  turn. 

The  final  turn  angle  has  been  given  by  Equation  (21)  as 


(B-8) 


Here  we  drop  the  subscript  0  because  Equation  (B-8)  can  now  be  applied  to  a 
spherical  earth  also.  In  the  main  text,  most  computations  were  made  for 
t!f£  <  90°  .  When  C  is  very  small,  however,  the  final  turn  angle  can  be  greater 
than  90  degrees.  If  ^  =  180°  ,  the  vehicle  would  have  made  a  lateral  displace¬ 
ment  of  2C  with  no  forward  displacement  and  would  fly  from  that  instant  on  in 
the  opposite  direction  to  the  original  heading.  If  il/^  =  360°  ,  the  vehicle  would 
have  made  a  full  360-degree  turn  and  find  itself  back  to  the  original  location  at 
a  much  lower  velocity.  In  such  case,  the  vehicle  would  be  dissipating  its 
kinetic  energy  through  drag  while  turning  and  it  would  have  no  loss  in  altitude 
and  no  gain  in  either  the  longitudinal  or  lateral  displacement  at  the  end  of  the 
turn. 

While  maneuvers  of  180-degree  and  360-degree  turn  appear  to  be  extremely 
interesting,  the  obvious  question  that  has  to  be  answered  is  whether  or  not  it  is 
practical.  To  answer  this  question,  we  examine  a  vehicle  with  L/D  =  2  , 

W / C^A  =  25  lb/ ft^  ,  at  V.  =  32,  000  ft/  sec  .  In  order  to  have  =  180°  , 
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the  turn  must  be  performed  at  an  altitude  of  184,  000  feet.  For  a  360-degree 
turn,  the  altitude  must  be  104,000  feet.  The  corresponding  value  of  C  is 
0.  0704  and  0.  00282,  respectively.  The  radius  of  the  turn  becomes  242 
nautical  miles  and  9.  7  nautical  miles,  respectively.  The  centrifugal  force 
created  by  these  turns  is  in  the  tens  of  gfs  for  ^  =  180°  and  in  the  hundreds 
of  g's  for  =  360°  .  For  lower  L/D  ,  higher  W/C^A  and  lower  V\  ,  it 
becomes  more  difficult  to  perform  such  turns  in  the  sense  that  the  altitude  will 
be  lower  and  the  centrifugal  force  will  be  higher.  Thus,  one  is  led  to  conclude 
at  this  time  that  constant  altitude  turns  at  very  small  values  of  C  are  not  likely 
to  be  practical. 

The  above  discussion  does  clarify  the  statement  made  in  the  text  about 
the  incorrectness  of  the  curves  of  lateral  range  and  down  range  at  very  sm.all 
W / (Cj^A  p)  as  shown  in  Figures  5  and  8.  As  one  can  see  from  Figure  B-l, 
the  lateral  range  and  down  range  must  have  the  characteristics  of  a  damped 
oscillation  as  W/(C^A  p)  goes  to  zero.  This  is  schematically  shown  in 
Figure  B-2.  Because  of  the  smallness  of  the  region  in  which  the  damped 
oscillation  character  occurs,  it  was  not  shown  in  Figures  5  and  8  to  avoid 
confusion. 


Page  48 


\f/  =  0°,  360°  - •- 

I  ORIGINAL  HEADING 

OF  VEHICLE 


Figure  B-l.  Vehicle  Flight  Path  at  Very  Small  C  . 


C0A  p 


Figure  B-2.  Schematics  of  the  Variation  of  \  and  (jl 
with  W/(CqA  p)  at  Very  Small  C  . 


NOMENCLATURE 


A  s  reference  area 

B  =  abbreviation  for  W/(C^A  p) 

C  =  abbreviation  for  2  B/(L/D) 

CD  =  drag  coefficient 

D  =  drag  force 

g  =  gravitational  acceleration 

L  =  aerodynamic  force  perpendicular  to  drag 

t  =  time 

V  =  velocity 

Vs  =  orbital  velocity  (=  26,000  ft/sec  for  earth) 

W  =  vehicle  weight 

X  =  abbreviation  for  In  Y 

\  =  lateral  range  angle 

|i  =  down  range  angle 

p  =  atmospheric  density 

0  =  bank  angle 

t|r  =  turn  angle 

Subscripts 

0  =  "flat  earth’1  solution 

1  =  spherical  earth  correction 

f  s  at  the  end  of  the  constant  altitude  glide 
i  =  initial  condition 

Superscripts 
* 


=  maximum 
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